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Abstract 

In this paper, we present the geodesic-like algorithm for the computation of the shortest path 
between two objects on NURBS surfaces and periodic surfaces. This method can improve 
the distance problem not only on surfaces but in R 3 . Moreover, the geodesic-like algorithm 
also provides an efficient approach to simulate the minimal geodesic between two holes on 
a NURBS surfaces. 
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1. Introduction 

Computing the distance between two objects on a surface plays an important role in 
many fields such as CAD, CAGD, robotics and computer graphics etc. In the Euclidean 
3-space M 3 , it has a simple mathematical presentation as following: 

min ||p — q\\ , (1) 

pSci, q£C2 

where C\ and c<± are two objects in M 3 . 

Although this representation is simple, however, it is hard to improve in general. The 
simplest case is the distance between two points and it can be estimated exactly by the 
Pythagorean theorem. If only one object is a point, then this problem is equivalent to the 



orthogonal projection problem, which has many applications [18L |19|. Many investigators 



have investigated the orthogonal projection problem. Chen et al. [3|,|4|, MaYL et al. [15| and 



Selimovic et al. [25| presented some effective methods that improve the distance problem 
between a point and a NURBS curve. Hu et al. jl] developed a good method to improve the 
orthogonal projection onto curves and surfaces. For the case that none of these objects is a 



single point, Kim [12[ presented a method to estimate the distance between a canal surface 



* Corresponding author: S.-G. Chen 
Email addresses: whclien@triu.edu.tw (Wen-Haw Chen), csg@mail.ncyu.edu.tw (Sheng-Gwo Chen) 



Preprint submitted to Elsevier 



November 25, 2010 



and a simple surface in 2003, while Chen et al. [5| improved this problem on two implicit 
algebraic surfaces in 2006. 

Unlike one can find many methods to investigate the distance problem in M 3 , there are few 
methods to study that on a curved surface. Maekawa [16j presented a very good method for 
solving the shortest path and the orthogonal projection problems on free-form parametric 
surfaces. Generally, the distance problem on a regular surface is more complicated than 
that in M 3 , even though the distance is just between two points on the surface, which 
is equivalent to find the length of the shortest path between them. We can find more 



information in reference [17|. This classical problem has many applications, such as in 



object segmentation, multi-scale image analysis and CAD etc. |2|. Il3l l24j . There are also 



■Y_' 
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many methods to estimate the shortest path on triangular mesh [14], [26j, polyhedral |10l 121 
and regular surface [HI, [23| etc. These methods can be extend to improve the distance 
between one point and one curve or between two curves on surface but they are not effective 
methods. 

In 2009, Chen j^] presents a new method to find geodesies on surfaces by the system of 
geodesic-like equations 

VE(u llVj ) = 0. (2) 

In this paper, the distance problem on NURBS surfaces and parametric surfaces will be 
improved by the geodesic-like method with B-spline basis. In fact, the geodesic-like method 
can also estimate the distance between two objects in M 3 but its efficiency is less than other 
algorithms that we known. 

This paper is structured as follows. Section 2 describes the definition of the distance 
problem on regular surfaces and the notion of geodesic-like curves. We shall present our 
geodesic-like algorithm to estimate the distance between two objects on surfaces, especially 
on periodic surfaces, in section 3. Section 4 presents some examples about the distance 
problem on NURBS surfaces by simulations. Finally, we illustrate a discussion about our 
method and conclude this paper in Section 5. 

2. Preliminaries 

Let us introduce the distance problem between two curves on a regular surface and the 
system of geodesic-like equations in this section. Suppose that S is a regular surface and 
(U, x) is a system of coordinates on S. A curve j(t) = (xj (t), X2(t)) is a geodesic curve in 
(U, x) on S if it satisfies the system of geodesic equations [l( 

From calculus of variation, a geodesic has an equivalent definition as below. 

Definition 1. A geodesic on a regular surface is a critical point of the energy variations. 
That is, the geodesic j(t) is a critical point of the energy function 

E(s) = lf\\%(s,t)\\ 2 dt, (4) 
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where f(s,t) is any proper variation of 'j(t). 



Here is a basic relationship between the length and energy functions. 

Theorem 2. Let S be a regular surface andp, q G S be two distinct points. If a is a shortest 
path between p and q on S, then a is a geodesic on S which pass through p and q. That is, 
the geodesic j(t) is a critical point of the length function 

L{s) = j\\?l{s,t)\\dt. (5) 

The distance between two points on a surface S is defined by the length of minimum 
path on S from p to q. Then 

d{p,q) = minL(7), (6) 

where V is the set of all paths on S from p to q and L{^) is the length of the curve 7 on S. 
From Theorem |2l the set T can be only considered the set of all geodesic on S from p to q. 

Consider the parametric surface S with a parametrization x : U — > M 3 and two curves 
Ci and c 2 on S. For simplicity, we denote Ci,c 2 : [a, b] — > U such that ci = x(ci([a, &])) and 
c 2 = x(c 2 [a, b])). Thus the distance between c\ and c 2 on S can be computed by 

d(ci,c 2 )= min d(x(ci(s)), x(c 2 (t))). (7) 

s,te[a,b] 

That is, d(ci, c 2 ) is the length of minimal geodesic from ci to c 2 . Equation (J2J) introduces 
a simple algorithm to improve this distance problem but it is too expansive. Let us describe 
it roughly. 

Algorithm 1. First, we digitize the curves Ci and c 2 to two sequences of points, {p«}^i 
and {qj}^ = Q, respectively. For each i,j, estimating the minimal geodesic 7^ between p$ and 

qj. Then the shortest path in { 7ij } (7^)^ (0 0) a PP r oaches the minimal geodesic between ci and 
c 2 on surface S when m, n are large enough. Of course its length approaches the minimum 
distance between Ci and c 2 on S. 

Solve the geodesic between two fixed points is crucial to solve Algorithm [TJ One can find 
many effective methods in the references js), HH, 14, 21, 23, 26]. However, if the numbers 



of {p^ and {qj} are large, this algorithm becomes very slow. In fact, Algorithm [T] is the 
simplest and the slowest method to improve this problem. 

We will improve the geodesic problem by the notion of geodesic-like curves jsf. 

Definition 3. Let x(w, v) be a parametrization of a regular surface S , x : U C M 2 — )■ S . A 

curve a(s) on U is called a geodesic-like curve of order n+l on S ifa(s) = ^™ =0 A^ n (s)(Ctj, Vi) 
is a B-spline curve and satisfies the system of geodesic equations 

(VE)(u l ,v J ) = 0, (8) 
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where 



1 r b 



E(u l ,v J ) = - j \\x(a(t))\\ 2 dt 
is the energy function of curve 

n 



8=0 



and (VE)(ui, Vj) is the gradient of E(ui,Vj) . 

Equation ([8]) is called the system of standard geodesic-like equations. Although the sys- 
tem of geodesic-like equations are integral equations, they can be improved by the Newton's 
method, the iterator method or other numerical methods 0, EH, 20, 22, 27] effectively. 



Since any piecewise differential curve can be approximated by the B-spline curves, a 
geodesic-like curve approaches a geodesic on S when the order of the geodesic-like curve is 
large enough. In the other words, we can estimate the distance between two points on S via 
the minimal geodesic-like curves. We summarize this property as follows. 

Theorem 4. Let S be a parametric surface and let 7 : [0, 1] — > S be a geodesic. Assume 
that the curve a n = 'YM =0 N™{t)(ui 1 Vi) is the geodesic-like curve between 7(0) and 7(1) for 
each positive integer n > 2. Then 

lim a n = 7. (9) 



3. Distance problem by geodesic-like algorithm 

The system of geodesic-like equations provides an elegant method to improve the distance 
problem between two objects on surfaces. We are now in a position to introduce this method 
in this section. The parametrization x on S is defined on U = [a, b] x [c, d\. That is 

x : [a, b] x [c,d] ->• S C M 3 . 

Let Ci and C2 be two differentiate parameterized curves on S and 

ci(s) : [0, 1] -> [a, b] x [c, d] 
c 2 (t) : [0,1] ->• [a,b] x [c,d]. 

Thus Ci = x(ci([0,l])) and c 2 = x(c2([0,l])) are two curves on S. To exclude the zero 
distance case from our consideration, we can assume that the two curves have no intersection. 
Denote Ci(s) = (u (s), v (s)) and c 2 {t) = (u n (t) , v n (t)) where u , u n : [0, 1] — > [a,b] and i> , 
v n : [0,1] — > [c, d] are all differentiable functions. Note that a B-spline curve a from [0, 1] 
into [a, b] x [c, d] with a(0) G C\ and a(l) G c 2 always has the form as 

«0) =Y^ZlN^x){u u v^ + N^x) Cl {s) + K{x) C2 {t) 

= J27^Nr{x){u i ,v i ) + N-(x)(u (s),v Q (s)) +NZ(x)(u n (t),v n (t)) 1 ; 
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where x G [0, 1]. 

Hence, we rewrite the system of geodesic-like equations to the following three differ- 
ent forms. These formulas improve the distance between two curves on S, the orthogonal 
projection problem on S and the shortest path between two points on S, respectively. 

The system of geodesic-like equations between two curves: From the equation (fTO]) , 
the parameters of the energy function E are s, t, ui, • • ■ , u n -i,Vx, • • • , f n -i- The system 
of geodesic-like equations between two curves can be rewritten as 

(VE) = (E a , E t , E Ul , E U2 , • • - , E Un _i, E n ,E V2 , • - - , £^„-i) = (11) 

The system of geodesic-like equations between one point and one curve: If c x is 

a constat curve on S, then the derivative of E about t is vanish. Thus we obtain the 
geodesic-like equation between one point and one curve. 

(V.E) = (E t , E U1 ,E U2 , ■ ■ ■ , Eu n _i, E V1 , E V2 , • • • , E Vn ^{) = 0. (12) 

Of course, The orthogonal projection projection problem on surface cab be improve 
by equation ffT2"]) . 

The system of geodesic-like equations between two points: Moreover, if ci and c 2 
are both constant curves on S, then the geodesic-like equations between these two 
points is 

(VE) = (E Ul ,E U2 , ■ ■ ■ , E Un _i, E Vl , E V2 , ■ ■ ■ , E Vn -i) = 0. (13) 

A curve satisfies one of equations ffTTl) - ffTSl is called a geodesic-like curve between ci 
and c 2 . Let us describe how to find the local minimal geodesic-like curve between two curves 
Ci and c 2 on the surface S. In this algorithm, we solve the system of geodesic-like curve 
equations by the Newton's method and the iterator method. 

Algorithm 2. (Geodesic-like algorithm) 

Step 1: Given two closed curves Ci and c 2 on the surface. Input an initial curve a such 
that the endpoints of a are on Ci and c 2 . 

Step 2: Solving the geodesic-like equations (equation / TJ7]) or / T73j) ) by the initial curve a 
and obtain a geodesic-like curve, which we still denote it by a, between c\ and c 2 . 

Step 3: If the set (a D Ci) U (a PI c 2 ) consists of the endpoints of a, then a is the local 
minimal geodesic-like curve between Ci and c 2 . Otherwise, trimming away some parts 
of the curve a such that the intersections of this trimmed curve, which we still denote 
it by a, and (ckDcx) U (aflCa) are only the endpoints of this trimmed curve (see Figure 
C]). Then repeat Step 2. 

By Theorem HI one will proceed by the geodesic-like algorithm to obtain the shortest 

path between ci and c 2 when the order n is large enough. We summarize it as follows. 
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Trimming 



Figure 1: the method of trimmed curves 

Theorem 5. Let S be a parametric surface and c 1; c 2 be two closed curves on S . For each 
n >2, a n is the local minimal geodesic-like curve that obtained by the geodesic-like Algorithm 
(algorithm^). If the set {a n } is a convergent sequence, then there exists a local minimal 
geodesic 7 between Ci and c 2 such that 

lim a n = 7. (14) 

n— ¥00 

Moreover, orthogonal to Ci and c 2 when n is large enough. 

3.1. Periodic surfaces 

If S is a periodic surface about one or two directions, then we may not obtain the 
local minimal geodesic-like curve by the Algorithm [2j It is because that the minimality 
of geodesies may not be preserved by the map of parametrization. To avoid this problem, 
we shall rewrite the domain U of parametrization of S. For simplicity, we assume that 
the original domain of parametrization is U — [0, 1] x [0, 1] and then the parametrization 
x : [0, 1] x [0, 1] ->• S. 

One-directional periodic surface Assume that S is a parametrization surface with u- 
directional period and x is the parametrization on S. Its domain of parametrization, 
the -u-u-plane, is as in Figure 2. Then x(0, v) = x(l, v) for each v G [0, 1]. Hence there 
is a function x : R x [0, 1] — > S such that x(tt, v) = x(m, v) for some u G [0, 1] provided 
(u — u) is an integer. Using the map x, we can find two curves c\ and c\ from [0, 1] to 
R x [0, 1] such that x(c?([0, 1])) = x(cj([0, 1])) = ci on S. Moreover, we assume that 
c\ = ci C [0,1] x [0,1] and c\ C [1,2] x [0,1]. Using Algorithm HI the geodesic-likes 
curves from c° to c 2 and from c\ to c 2 on U can be found and we denote them by 70 
and 71 respectively (see Figure [2]). Then the one in {x(7o), x(7i)} with smaller length 
is the local minimal geodesic-like curve between ci and c 2 on S. 

Two-directional periodic surface If S is a periodic surface about two directions, then 
we shall find x : R x R — > S such that x(u, v) = x(w, v) for some u, v G [0, 1] if both 
{u — u) and (v — v) are integers. Similarly, we can find four curves C [0, 1] x [0, 1], 
c\ C [1,2] x [0,1], cj C [0,1] x [1,2] and c\ C [1,2] x [1,2] (see Figure ED such that 
x(c^) = x(ci) = Ci on S for i = 0,1,2,3. Denote the minimal geodesic-like curve 
between c\ and c 2 by 7$ for i = 0, 1,2,3. Thus the minimal geodesic-like curve on S 
between ci and c 2 on S is the one in {x(7i)}f =0 with the shortest length. 
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Figure 2: Local minimal geodesic-like curve on a one-directional periodic surface. 
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Figure 3: Local minimal geodesic-like curve on a two-directional periodic surface. 
4. Simulations 

To apply our method in practice, we present some examples by simulation. The geodesic- 
like curves in our simulations are all uniform quadratic B-spline curves in R 2 . 

First we consider an open surface S and two closed curves c x and c 2 on S as in Figure HI 
The surface S is a cubic B-spline surface with (8,4) control points. The red curve in figure 
4 is the local minimal geodesic-like curve of order 11 between Ci and c 2 and its error is less 
than 10" 6 . 

Secondly, a surface of revolution is an example of one-dimensional periodic surfaces. 
Figure is the domain of parametrization (wf-plane) of the surface of revolution as in 
Figure El In Figure El there are two geodesic-like curves in the wv-plane, one is from c\ to 
C2 and the other is from c\ to c 2 . Then the image under parametrization of the shorter one 
is the local minimal geodesic-like curve between these two curves on the surface. 

Thirdly, a typical example of two-dimensional periodic surfaces is the torus. Figure [7] is 
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Figure 4: Distance between two closed curves on a NURBS surface 
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Figure 5: domain of parametrization of a u-periodic periodic surface 



the domain of parametrization of a torus as in Figure [SJ There are four geodesic-like curves 
in the ttf-plane. Therefore, the image under parametrization of the shortest one is the local 
minimal geodesic-like curve between two closed curves on the torus. 

Lastly, we construct a face model as in Figure [9] by NURBS surface and find the minimal 
geodesic-like curves between two holds (the eyes) on the surface. The data in FigurdTUl are 
about the geodesic-like curves of different orders between the two holes in Figure 9. Here in 
Figure [TD] the order means the number of control points while error (%) is the percentage of 
error, which is defined by 



error (%) 



Length — minimum distance 
minimum distance 



x 100%. 



(15) 



The red curve in Figure is the local minimal geodesic-like curve of order 30 and the green 
curve is the exact minimal geodesic between two holes. Then the lengths of geodesic-like 
curves constructed by our method approaches the minimum distance between the two holes. 
To deserve to be mentioned, error(%) will be less than 10~ 7 provided the geodesic-like 
curve is constructed by 60 control points. It proposes that the geodesic-like algorithm has 
increased actually computational efficiency of this simulation . 
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Figure 6: One-directional periodic surface 
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Figure 7: Distance between two closed curves on torus (uv-plane) 

5. Discussion 

The geodesic-like algorithm provides an effective and reliable computation of shortest 
paths between two curves on surfaces. For computing the shortest paths between two curves 
on R 3 , our method is comparable with other well-known methods. Especially, the construc- 
tion of geodesic-like curves only bases on the uniform quadratic B-spline curves since it is 
enough to us to consider the geodesic-like curves in the plane. Significatively, our method 
can be extended to solve the distance problem between any two objects on surfaces and the 
distance problem in higher dimension. 

To solve the system of geodesic-like equations, however, Newton's method is too expan- 
sive. Moreover, it can only solve local minimal geodesic-like curves but not global minimum 
ones. In the future investigation, we expect to find a numerical method to solve efficiently 
all local minimal geodesic-like curves between two objects on surfaces to overcome these 
problems. 
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Figure 8: Distance between two closed curves on torus 




Figure 9: The distance between two holds on a face model 
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